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Abstract 

We use the Bogoliubov formalism to study both, particles and gravitons creation at the reheating 
epoch, after a phase transition from inflation to a radiation dominated universe. The modes of 
the inflaton field fluctuations and the scalar fluctuations of the metric at the end of inflation are 
obtained by using a recently introduced formalism related to the Induced Matter theory of gravity 
The interesting result is that the number of created particles is bigger than 10 90 on cosmological 
scales. Furthermore, the number of gravitons are nearly 10 -17 times smaller than the number of 
created particles. In both cases, these numbers rapidly increase on cosmological scales. 
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I. INTRODUCTION 



Inflation has become the standard paradigm for explaining the homogeneity and the 
isotropy of our observed Universe [l, 2|. During this epoch the energy density of the Uni- 
verse was dominated by some scalar field (the inflaton), with negligible kinetic energy density, 
such that the corresponding vacuum energy density was the responsible for the exponential 
growth of the scale factor of the universe. During this phase a small and smooth region of 
the order of size of the Hubble radius grew so large that it easily encompassed the comoving 
volume of the entire presently observed Universe. This is the reason for which the observed 
Universe is so homogeneous and isotropic. Furthermore, it is now clear that structure in 
the Universe comes primarily from an almost scale-invariant super-horizon curvature pertur- 
bation. This perturbation originates presumably from the vacuum fluctuation, during the 
almost-exponential inflation, of some field with mass much less than the Hubble parameter 
H . Indeed, any scalar field whose mass is lighter than H suffered, in a (quasi) de Sitter 



epoch, with a scale independent quantum fluctuations spectrum 

During inflation, particle production can only occur for particles that are light compared 
to the Hubble scale without classical conformal invariance; gravitons and massless minimally 
coupled scalars and light fermions are unique in that respect. Particle creation during 
inflation was studied many years ago in the framework of standard inflation {7, 8]. The process 
appeared to be straightforward: in models like new inflation and chaotic inflation 9J. These 
models incorporate a second order phase transition to end inflation, the inflaton field would 
wind up oscillating around the minimum of its potential near the end of inflation. These 
oscillations would produce a sea of relativistic particles, if one added (by hand) interaction 
terms between the inflaton and these lighter species. The direct production from vacuum 



fluctuations during inflation of X-particle was considered in 10 1. Particle creation has been 



also considered in warm inflationary 



ll| and fresh inflationary 



121 ] scenarios. More recently, 



particle creation during inflation was considered in a model of brane inflation where two 
stacks of mobile branes are moving ultra relativistically in a warped throat 13 ]. 

After inflation, the universe could have suffered a phase transition from a de Sitter vacuum 
dominated state to a decelerated radiation dominated stage. At this moment, a great amount 
of the potential density energy is transferred to radiation energy density due to the decay of 
the inflaton field produced by its interaction with other boson fields 14Hl8|. Some years ago 
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Kaiser[19| demonstrated that particles produced from the parametric resonance effect when 
H ^ is fewer than in the H ~ case, but can still be exponentially greater than when the 
resonance is neglected altogether. 

On the other hand, in a previous work [2^] we studied the scalar metric fluctuations of a 5D 
spacetime background metric (Ai,g), which is Riemann-flat and hence Ricci-flat: Hab — 0. 



From the mathematical point of view, the Campbell-Magaard theorem [2 lH24} serves as a 
ladder to go between manifolds whose dimensionality differs by one. This theorem, which is 
valid in any number of dimensions, implies that every solution of the 4D Einstein equations 
with arbitrary energy momentum tensor can be embedded, at least locally, in a solution of 
the 5D Einstein field equations in vacuum. Because of this, the stress-energy may be a 4D 
manifestation of the embedding geometry. Physically, the background metric there employed 
describes a 5D extension of an usual de Sitter spacetime. Inflationary cosmology can be 
recovered from a 5D vacuum 
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27(| . Other version of 5D gravity, which is mathematically 



similar, is the membrane theory. In this theory gravity propagates freely on the 5D bulk 



and the interactions of particles are confined to a 4D hypersurface called "brane" 28M30 1 . 
Both versions of 5D general relativity are in agreement with observations. 

In this work we study both, particle and graviton production after a phase transition 
occurred after inflation. To make it, we shall consider the solutions obtained at the end of 
inflation from a 5D vacuum on an effective 4D hypersurface described by a de Sitter metric. 
After it, we shall use the Bogoliubov formalism to calculate the modes solution for the fields 
after a phase transition to a radiation dominated epoch. The topic here studied is very 
important, because the case of inflaton decay into further inflatons may be of interest for 
dark matter searches. Following their production, the inflatons would decouple from the 



rest of matter 



31] . If the inflatons were given a tiny mass, then these bosons could serve as 



a natural candidate for the missing dark matter 19]. 



II. INDUCED MATTER AND EMBEDDINGS 



We consider a 5D manifold (Ai, g a b) with a coordinate system x = {x a }. We are interested 
in a 5D theory of gravity on which we define a 5D vacuum, such that the first variation of 
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the action is ^5X =^ 5X E +^ 5X M , where 

The first term is the variation of the gravitational Einstein action ^X E and the second one 
is the variation of the matter action ^Im- Here, G is the gravitational constant, g is the 
determinant of the covariant tensor metric g ab l and R is the 5D Ricci scalar on the metric. 
The energy-momentum tensor of matter, T ab is defined from the variation of the matter 
action ^Xm under a change of the metric, and will be considered as null to describe the 5D 
apparent vacuum. The Einstein equations on the 5D manifold M. (we use c — h—1 units) 

G\ = -MGfl (2) 

where G% and T b a are respectively the Einstein and Energy Momentum tensors on M.. Fur- 
thermore, we can define a scalar function l(x), which represents the foliation of the higher- 
dimensional manifold. We shall consider that the extra coordinate is space-like. Each 
hypersurface, Sj, is considered as a 4D Lorentzian spacetime. We denote n a as the normal 
vector to the hypersurface E; 2 

n a = -VJ, (3) 

such that the expression n a n a = — 1 normalizes n a . Now we can define a coordinate system 
y = {y a } on £/. The basis vectors are 

3x a 

< - 5?. - 0- ( 4 ) 

These objects can be used to project 5D tensors (for instance, the metric tensor) into 4D 
ones (which lives on the Ej hypersurface) 

e a a e b p g ab = e^ej (h ab + n a n b ) = h a p. (5) 



1 In this work a, b run from to 4 and Greek letters rum from to 3. 

2 In what follows we shall denote the covariant derivative on the 5D hypersurface as V a (•■) and the covariant 
derivative on with a semicolon: (..) ;o . 
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A. Einstein tensor on Sj 

The extrinsic curvature K a p of the 4D hypersurface E/ is a symmetric 2-range tensor 
given by the derivative of the induced metric in the normal direction to E ; 

K aP = e a a e b p V a n b . (6) 

If now we consider an alternative coordinate system on M. : {y a , I}, such that 

dx a = e a a dy a + l a dl, 

where the vector tangent to lines with constant y a can be decomposed into the sum of a 
part tangent to Ej, and a part normal to Ej 

T = N a e a a + n a , 

such that l a d a l = 1. The 4D vector N a is the shift vector, which describes how the y a 
coordinate system changes as we move from a given Ej hypersurface to another. Finally, the 
5D line element dS 2 can be written as 

dS 2 = h ab dx a dx b = Kn [dy a + N a dl) {dy p + N?dl) - dl 2 , (7) 

which, for a constant foliation dl = 0, reduces to 

ds 2 = h aP dy a dyP. (8) 

In this work we are interested in dealing with a 5D Riemann-flat metric, so that the Ricci 
tensor is null 

Rob = 0. (9) 
On each E; hypersurface, the Gauss-Codazzi equations are 3 

Rabcd e a e /? e 7 e a = ~ 2K a [sK^]p, (10) 

R mabc n m e a a ey i = 2K a[ p. n] . (11) 
If we use the expression for the 5D Ricci tensor 

R ab = {h^e^el - n m n n ) R ambn , (12) 

3 For a 5D Riemann-flat metric, are Rap-yS = 2K a ^K^p, and K a ^.^ = 0. 
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we obtain the contractions of 



Rab C a 6g 


= 0, 


R ab e a a n b 


= 0, 


R ab n a n b 


= 0. 



(13) 



Using ffl2l) into f|T3|) . and making use of ffTTl) . we obtain the expressions 



^ = E aP + K£ (Kfr - K h Pll ) , (14) 
(p<*P) a = [K afi - h ap K) , a = 0, (15) 
E^hT = 0, (16) 

where K = h a ^K a p, E a p = Ep a and E a/3 = R manh n m e a a n n Notice that the equation (115]) 
means that the second rank (symmetric) tensor P a " is conserved on the 4D hypersurface 
£j. It is important to notice that the Ricci tensor on related to the 5D Riemann flat 
metric R mn io = 0, is given by 

R afj = K£(K f)li -Khf ) r). (17) 
Finally, the Einstein tensor G Q/3 = R a/3 — g a/3 R/2 on a given £j, is 

= + R «prf _ fffi K l* ( lg ) 

Notice that for a 4D metric induced from a 5D Riemann-flat metric, the Einstein tensor 
only depends on the extrinsic curvature 

G afi = K$prt - -h^KT^. (19) 
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B. Einstein equations on 

In this work, we are interested in describing the primordial universe on a 4D hypersurface 
Sj. In order to make it, we shall consider a massless single scalar field f{x a ), which, we shall 
consider to describe the physical vacuum on the 5D manifold Ai. The kinetic Lagrangian 
corresponding to this field is 

g^aVfi, (20) 
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where g = I 8 exp 6(iV — 4ip) is the determinant of the covariant metric tensor and go is a 
dimensional constant. Hence, to describe the energy momentum tensor related to (p, on A4, 
we set 



-8nGf ab - 9Cip 







dg, 



ab 



dx r 



d(g, 



ab,m ) 



(21) 



It is easy to demonstrate that the T ab projected on the hypersurface £/ is 



rpa/3 rpab a 8 

1 — 1 ^a^b- 



(22) 



Finally, using ( 122]) and (fl9|) . we obtain the Einstein equations on obtained from a 
Riemann-flat 5D vacuum 



G af}\ 



f-l 9 fw 



1 dC^ 



e a e b 



d 



dg, 



ab 



dx T - 



_d(g a b,r 



(23) 



that describes the equations of motion for the fields on hypersurfaces with constant I (i.e., 
for constant foliations). Notice that the equations of motion obtained from a Ricci-flat 5D 
metric must include the Einstein tensor (fl8|) , rather than (Tl9|) , in the left side of the Einstein 
equations fl23|) . 



III. REVISITING THE FORMALISM FOR SCALAR METRIC FLUCTUATIONS 
ON A 5D RIEMANN-FLAT METRIC 



We consider the Riemann-flat background metric 32| 



dSl = l 2 dN 2 - l 2 e 2N dr 2 - df 



(24) 



where dr 2 = 5ijdx l dx^ . Here, x % are the 3D cartesian space-like dimensionless coordinates, iV 
is a dimensionless time-like coordinate and / is the space-like non-compact extra coordinate, 



which has length units. The non-perturbative metric fluctuations of t 



(|24|) . are introduced in our analysis by the line element introduced in 20] 

dS 2 = l 2 e 2 ^dN 2 - l 2 e 2(N -^dr 2 - dl 2 , 



re background metric 



(25) 



where the metric function ij)(N,r,l) describes the gauge-invariant metric fluctuations with 
respect to the background metric. In order to describe a 5D physical vacuum, we shall 
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consider a massless and free test scalar field ip = <p(x a , I) defined on (|25|) . The dynamics of 
(p can be derived from the action 



d 4 x dl\ 
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R 



16nG 



(26) 



Now we consider a semiclassical approximation for the 5D scalar field ip in the form 
ip(N,r,l) = <pb(N,l) + 5<p(N,r,l), with (pb denoting the background part of ip and S(p 
denoting the quantum fluctuations of ip, which will be considered as very small. This is 
consistent with a linear approximation for the scalar metric fluctuations. Thus, a first-order 
approximation in the gauge scalar fluctuations of the form e ±2 ^ ~ 1 ± 2tp, will be sufficient 
in order to have a good description of these fluctuations during inflation, in the present 
formalism. After making this approximation, one find that the dynamics on 5D of (fb{N, I) 
and Sip(N, r, I), are 



d 2 ip b 3 d(p b 



8N 2 
d 2 5ip 



8N 
d5ip 



'9 2( fb , AJ d(p 



di 2 

2Nt-j2 



+ 41- 



81 



0. 



(27) 



- 3^-e^V z r 5ip 

\ 2 d\pb + Al 8ip h 



8N 2 ' " 8N 
- 2^ 



,8 2 5ip 



o. 



+ 41 



85ip 

~dT 



(28) 



dl 2 81 

The expression (|27|) gives the dynamics on the background of <p(N,r,l), whereas Eq. ([25]) 
describes the dynamics for the quantum fluctuations 5<p(N, r, I) in terms of the scalar metric 
fluctuations ip and the background field ipb- On the other hand, the expression for the energy 
momentum tensor components on the background is 4 



87rGT f 



8C lr 



8 



ab 



8g 



ab 



8x r 



8C„ 



d (g ab m ) 



(29) 



9ab=9ab, l P= l Pb 

where the Lagrangian density corresponding to the inflaton field is 
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g ab y, a y,b- 



(30) 



One can find the linearized equation which describes the dynamics of the salar metric 
fluctuations ip on the linearized fluctuating metric dS 2 = (1 + 2tp)l 2 dN 2 — l 2 e 2N (1 — 2tp)dr 2 — 



We denote with g a b the background tensor metric and ifb the solution for ip(x a ) on g a b. 
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# 2 |20|: 



14Z^ + 4f 
dl 



,d 2 *p 



dN dN 



d(fi 
dl 



2 d(pb d8y 
+ dl dl 



where ip complies with the linearized non-diagonal Einstein equations 



d 2 ip dtp 



+ 



dx i dN dx 



AttG 



difb d5ip 
dN dx i ' 



d 2 ip dtp 8nG ( d5ip dip b d(p b d5(p 
dldN + ~dl = ~Y~ \dN dl + dN dl 

d 2 ip _ ^ G ( difb_d5ip_ d5(p d(p b 
dx % dl V dx % dl dx % dl 



(31) 

(32) 
(33) 

(34) 



These equations provide us with the dynamics for both, the scalar field and the 5D gauge- 
invariant metric fluctuations}^ . 



IV. DYNAMICS OF FIELDS AT THE END OF INFLATION ON AN EFFECTIVE 
4D DE SITTER BACKGROUND 



We consider the background metric (124"]) . If we take a static foliation 1 = 1$ = 1/Hq, with 
the transformations 



N = t/l = H t, 
R = l r = 1/H r, 

we obtain the effective 4D de Sitter background metric 

dS 2 = dt 2 - e 2Hot dR 2 . 



(35) 
(36) 

(37) 



Here, H Q is a constant that represents the Hubble parameter. For the relativistic point of 
view, this means that now we shall move with penta-velocities 5 

L, u l = u l = 0. 



u 



The effective 4D action on £# derived from the 5D action ( 126]) . reads 



(4) T 



eff 



d x 



( 4 )< 



(4) 



90 



V 16ttG 2 



+ -g^,^, v + V 



(38) 



(39) 



Latin indices denote 3D spatial coordinates. 
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where $(t, r) = <p(t,r,l)\ l=H -i is the massive scalar field induced on E#. Furthermore, 
( 4 )g is the determinant of the 4D induced metric, which in the case of background metric 
([37)1 yields = exp(6H t). Moreover, in the case of the perturbed metric (125)1 gives 
( 4 )<7 = exp[2(3//o^ ~~ 2^)], and ^go is a dimensionless constant. The induced 4D effective 
potential V in the action (139 p . has the form 



V 



dip 



(40) 

N=H t,R=r/H ,l=l =l/H 

Given the quantum nature of the fields $ and \I> it seems convenient to use the quantization 
procedure. To do it, we impose the commutation relations 

$M),n° $) M')] = ^ (3) (R - R') , [*M),nj^ (*,#)] = ^ (3) (i? - i?') , (4i) 

where IT^ and IT^ are the canonical momentums for $ and respectively. We can derive 
the commutation relations on the effective 4D background metric. If the inflaton and metric 
fluctuations are small, these commutators can be approximated to 



$M),$Mo 



*M),*M') 




(42) 



(43) 



In this section we shall examine the evolution of the background inflaton field, jointly with 
the metric and inflaton fluctuations. 



A. Dynamics of the background field at the end of inflation 



The dynamics of the background inflaton dynamics on the effective 4D background (de 
Sitter) metric is described by the equation (127)) with I = l = 1/H and N = H t. In what 
follows we shall define = <fib(N = Hot, I — Iq — 1/H ). After making separation of 

variables in ( 12 7p . we obtain 



dt 2 



dt 



+ H 2 m 2 $ b = 0, 



m being a dimensionless separation constant. The background Friedmann equation G*. 
—8irG (T\), on the effective 4D hypersurface is 



d<fi 
dt 



H ^ l ~dT 



3 1 



l=H~ n 



AnGl 2 



3HI 
AnG' 



(45) 
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The complete solution for eqs. fl4*4|) and (l4"5j) is with m = 0, such that the inflaton field is a 
constant of time: $6 = W / ^^[20|. 

However, at the end of inflation the second term in the equation (jHj) can be neglected 
and that equation can be approximated to 

^ + Hlm^ b « 0. (46) 

In what follows we shall consider the dynamics of the inflaton field governed by the approx- 
imated equation ( l4"6l) with the effective 4D Friedman equation ( ]4"5i) . We shall consider that 
at the end of inflation the mass of the inflaton field is of the order of the Hubble parameter. 
In particular, if we consider m = 2, the dynamics of the metric fluctuations described by 
(|3T|) is simplified and we obtain the following solution for ( 146]) 

$ 6 (t) = $ cos(2H t + 6), (47) 

where <3>o = \\J~~^g anc ^ ^ * s some constant of integration to be determined. 

B. Dynamics of metric fluctuations at the end of inflation 

Using the result (14"?]) . and combining equations (T3TT) to (132]) over the 4D hypersurfaces 
Ejj : I = Iq = Hq 1 , we obtain an uncoupled equation for ty(t,r) 

— - + 3H — + 2H 2 ^ (1-3 cos 2 (2H t + 5) + S 2 tan(2H t + 5)) - 3H 2 e~ 2Hot V 2 r ^ = 0, 
ut at 

(48) 

where S 2 = 4kf, k\ being a separation constant. 

Following a canonical quantization process, the field $ can be expanded in Fourier modes 

3) = / Sk h e ^ (t) + fltfc e ~**s(*)] ' (49) 

where the annihilation and creation operators and a k , satisfy the usual commutation 
algebra 

[a k , al} = 6®(k - 9), [a h , a k ,} = [a\, a{] = 0. (50) 

Using the commutation relation ( 150]) and the Fourier expansion ( 149]) . we obtain the 
normalization condition for the modes 

mint) - zmut) = ^ (^) 3 ^ 3 , (5i) 
ii 



where a(t) = e Hot is the scale factor in the 4D metric and ao = e H ° to , being to the cosmic 
time at the end of inflation. Inserting (|49p in (j48j) . we obtain that the modes satisfy 



dt 2 



+ 3#o^ + 2# 2 £ fc ^1-3 cos 2 (2H t + 5) + ~k 2 e~ 2Hot + S 2 tan(2F i + <5) J = 0. (52) 



Lets us consider the phase 5 so that at the end of inflation 2H t + 5 ~ 7T. In this case, the 
equation (l52l) leads to 



<9 2 £ 



^2 " :u S • - , (53) 

By means of the equation floTj) and choosing the Bunch-Davies vacuum condition, the 
normalized solution of (1531) . is 

&(t) = ^VG(^) S/ V a) N(*)], (54) 
where 7i^[x(t)] is the second kind Hankel function, v = 5/2 and x(t) = y/3 k e~ Hot . 

C. Dynamics of the inflaton field fluctuations at the end of inflation 

Now, employing again the equations f l3Tj) to ( 152"]) over the 4D hypersurfaces : I = l = 
Hq 1 and incorporating the solution (JHJ), we obtain that the 4D quantum fluctuations for 
the inflaton field are determined by the equation 20] 



,d 2 5ip A1 d5(p 



dt 2 ' ~" u dt 
Performing a Fourier expansion of 5$ in the form 

1 



+ 41- 



di 



+ 8# 2 *$ fe = 0, (55) 



l=H~ n 



5<$>(t,R) 



(2tt) 3 / 2 



t „-%K R„ 



(56) 



where the annihilation and creation operators ax and a} K satisfy the commutation algebra 
given by floll . Inserting fl53|) and fT4T?j) in f[5"4|) . it leads to 



,00- %C , tt2 K 2-2H t tt2 



l 2 ^ + 4l dm 



01 



+ 8H 2 £ K $ b = 0. (57) 



The previous equation is inhomogeneous, non separable and presents a source term. The 
general solution will be the solution of the homogeneous equation plus a particular solution 
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VK (t) = [/(*) cos (kV3e- Hot ) + g(t) sin (W^T™)] + e — J— Uf [x(t)\ , (58) 

where f(t) and g(t) are polinomial functions of e~$~ , and [x(t)] is the second kind Hankel 



function with \i = ^ 9 2 4Af2 and x(t) = ke Hot . Here, we have used both, the Bunch - Davies 



2 

vacuum and the normalization conditions 



Vh* K VhK ~ VhKVhK = i (^) , ( 59 ) 

for the Fourier modes of the homogeneous solution. 

Towards the end of inflation one can neglect the lower terms of f(t) and g(t), to take 

H t 

only those with higher order in e 2 . The term with the Hankel function can be neglected 
too, so that we obtain 



'< t > = i ^ = d^V^ eiH °'° el,J "'' (60) 

Finally, the approximated solution for the modes of the inflaton field fluctuations at the end 
of inflation can be rewritten as 

VK(t) = f(t)e-^^ Ho \ (61) 

where f(t) is given by (1601) . 



V. PARTICLE CREATION AFTER A PHASE TRANSITION, AFTER INFLA- 
TION 

To consider the particle creation lets first rewrite the solution flBTl) in terms of the con- 
formal time in a spatially flat 4D Friedmann- Robertson- Walker cosmology dr = a~ 1 dt = 
e~ Hot dt. Furthermore, we obtain that the modes of the inflaton fluctuations are 

Vk(t) = f( T )e ikV ~ 3Hoao \ (62) 

and their derivatives with respect to r, which we denote with a prime, are 

rj' K {r) = f( T )e ikVEH °^ (zKV3H a - A) , (63) 
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where 



f(r) 



40 H el Hoto / y/3 



'-HoCLqt) 



(64) 



3/^(25 - 4/x 2 ) V 2K 5 

In order to calculate the number of particles produced for each mode (Nk), we match the 
solutions obtained at the end of inflation, with those of the long-wavelength modes during 
a phase transition from inflation to the radiation-dominated era. To calculate this last 



33j). For the long- wavelength 



solutions, we make use of the Bogolyubov transformation 
modes we can approximate the phase transition to be instantaneous, occurring at some time 
Tq. So, we have that 

for r < t : a(r) = — -4-, 

V ; Ht (65) 
for r > t : a(r) = — a (r — r), 

where r = tq — (H 2 tq)~ 1 . For convenience, we set o(to) = acb so we obtain that f = 0, 



Hn 



-1 



and Tq = —1. After the transition, the Hubble parameter is no longer constant, 



H = H(t), and we get: 



(66) 



for r < t : #(r) = H , 

for r > r : H(r) = ^. 
In what follows we shall compute the number of particles and gravitons in the radiation 
dominated era, which takes place after a phase transition from the inflationary epoch. 



A. Particle creation from inflaton field fluctuations in the radiation dominated 



era 



In order to calculated the number of created particles after a phase transition, we calculate 
the modes of the inflaton field fluctuations taking into account the Bogolyubov transforma- 
tions 



Vk{t > Tq) 



HqT 



so that their temporal derivatives are 



axe 



(67) 



V K i T > T o] 



Ho 



axe 



• (68) 
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To calculate the coefficients ax and (3k, we use the continuity of t]k and r)' K at r = To = — 1, 
so that 

(3 K = a /(r ) ^2KV3 (l - , (69) 

a K = 2 a /(r ) ^2^ - l) e" 2 ^, ( 70 ) 

where 



40 ffp et g °*° t V3 
/(r ° ) = 3M25-V)V2^- (71) 

Finally, by using Nr = \Pk\ 2 we obtain the number of created particles with a wavenum- 

ber K 

N K = 2a 2 |/(r ) | 2 (l + . (72) 

To estimate numerically the number of particles created after the phase transition on 
cosmological scales (small wavenumbers) , we shall consider that the number of e-folds, N e , 
at the end of inflation is at least 6 : N e > Hoto ~ 60. The value of K depends upon the 
details of the transition between inflation and the radiation dominated epoch. For long- 
wavelength modes (on cosmological super Hubble scales), the transition can be considered 
as abrupt. In general one expects that K < 1 34j. For instance, if we take 0.01 < K < 1 
and H ~ 6 x 10~ 6 M p 7 , with //(25 - 4/i 2 ) ~ 24 8 , we obtain 

5.5 x 10 78 <N K < 5.64 x 10 90 , (73) 

which, in cases where the wavenumber is small, is a number sufficiently important to agree 
with the desired values for entropy S > 10 89 . 

In the figure ([1]) was plotted the number of particles created for 0.01 < K < 1. Notice 
that the number of created particles Nk increases dramatically for cosmological scales. 



6 This means that we are considering that the phase transition occurs abruptly after inflation ends. 

7 The gravitational constant G = M~ 2 . 

8 This is the value corresponding to a scale invariant power spectrum (with /i = 3/2), for the inflaton field 
fluctuations on cosmological scales. 
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B. Graviton creation from scalar metric fluctuations 



We can also estimate the number of gravitons per mode generated by the metric fluc- 
tuations at the end of inflation. Again we first must rewrite the solution of the metric 
fluctuations modes in conformal time r. The modes for the metric fluctuations (1541) can be 
expressed 

6c(t) = C{T)e~ iK ^ Hoao \ (74) 
so that the conformal time derivative is 



&(r) = (C'(r) - iKV3H a ) e -**v^W 



where the function C(r) and its conformal time derivative are given by 

1 



C(t) 



1 /2;7G :e f^^a ll + 



1 



H$cl 2 K 2 t 2 



(75) 



(76) 



(77) 



3\/ y/ZK 

Let us consider again the phase transition (165]) . with a(r ) = a , so that we have f = 0, 
Ho = clq 1 and To = —1. After the phase transition, the metric fluctuations modes can be 
written as 



T 



A K e~ lK ^- 



B K e 



,KV3n 



so that 



A K e 



l + ^) + B K e^ 



(78) 



• (79) 



To obtain the coefficients Ak and Bk we use the continuity of £k and at moment of 
the phase transition, r = r = — 1, and we obtain 



3K 3 V 3 V 2 ' 



(80) 



B K = e 2«V3^ e fffoto; 



3K 2 



(K 2 - 2) + z f 1 - JL) 



V 3AT 2 / 



(81) 
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Using again Nfl = \Bk\ , we can obtain the number of gravitons generated after the 
phase transition 

N k = e*™*—^ _ _ + __j _ e^° V 3 V (K~* - 2K-*) 

+ ~Y~ (82) 

Notice that the third term inside the first bracket is the dominant one. Now, if we consider 
the values of the parameters at the end of inflation, we can calculate the approximate number 
of gravitons created after the phase transition. In order to estimate a numerical value, we 
set H t > 60 and K ~ 0.01, so that we obtain 

Njp > 2.5 x 10 73 , (83) 

which is a number very small with respect to the number of particles created (!73|) after the 
phase transition. In the figure ([2]) was plotted the number of gravitons created , which 
also increases dramatically on larger cosmological scales. 



VI. FINAL COMMENTS 



We have studied particle and graviton production after an abrupt phase transition from 
inflation to the radiation dominated era. We have considered the Bogoliubov formalism to 
calculate the modes solution for the fields after a phase transition to a radiation dominated 
epoch. Our approach is different from others developed in [l9|], because here resonance effects 
are not important for the inflaton and metric fluctuations. Furthermore, the solutions for the 
modes of the inflaton fluctuations and the scalar fluctuations of the metric were obtained 
using a recently developed formalism which is related to the Induced Matter theory of 
gravity. During inflation, such approach makes possible the simultaneous description of 
the coupled dynamics of the inflaton field fluctuations with the gauge-invariant fluctuations 
of the metric. After inflation this solutions can be matched with those of the radiation 
dominated epoch to obtain the Bogoliubov coefficients. The results here obtained are very 
interesting, because the number of created particles after inflation results to be greater than 
10 90 at cosmological scales, and increases dramatically as K~ 6 . Furthermore, Nk is very 
sensitive to the number of e-folds during inflation. This is due to the fact that iV^- increases 
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as a(t) 3 . In what respect to the number of gravitons created after inflation, the results are 
qualitatively similar to those of created particles. However, the number of gravitons created 
are, at least, of the order of 10 73 , which is much smaller (10~ 17 times) than N K . All these 
results were obtained using the values H t = 60, with H = 1CT 6 M p , and wavenumbers 
K < 10- 2 . 
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FIG. 2: Number of gravitons created after inflation , with respect to the wavenumber 0.01 < 
K < 1. 
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